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In a semimetal, both electrons and holes contribute to the density of states at the Fermi level. 

The small band overlaps and multiband effects engender novel electronic properties. We show that 
a moderate hydrostatic pressure effectively suppresses the band gap in the elemental semiconductor 
black phosphorus. An electronic topological transition takes place at approximately 1.2 GPa, above 
which black phosphorus evolves into a semimetal state that is characterized by a colossal positive 
magnetoresistance and a non-linear field dependence of Hall resistivity. The Shubnikov-de Haas 
oscillations detected in magnetic field reveal the complex Fermi surface topology of the semimetallic 
phase. Specially, we find a non-trivial Berry’s phase in one Fermi surface that emerges in the 
semimetal state, as evidence of a Dirac-like dispersion. The observed semimetallic behavior greatly 
enriches the material property of black phosphorus, and sets the stage for the exploration of novel 
electronic states in this material. 

PACS numbers: 71.18.-l-y, 71.30-l-h, 74.62.Fj 


The electronic band gap is an intrinsic feature of a 
semiconductor that largely determines the semiconduc¬ 
tor’s electronic and optical properties. In black phospho¬ 
rus, which recently emerged as a unique layered semi¬ 
conductor that exhibits high carrier mobility [H, 0, the 
puckered honeycomb lattice in each P atom layer [Sj dic¬ 
tates that the direct band gap is mainly determined by 
the out-of-plane p-,"like orbital of phosphorus and fur¬ 
ther reduced by the interlayer coupling in the bulk crys¬ 
tal 00. Thus, a slight change in the height of the 
puckered layers strongly modulates the direct band gap 
of black phosphorus Q. An ^plied hydrostatic pres¬ 
sure can close the band g^ [sj, III , and drive a structural 
transition above 4 GPa [8|. The effectiveness of such 
pressure-induced band modification and the theoretically 
suggested linear band dispersion along the armchair di¬ 
rection in black phosphorus 0 enable exploring exotic 
electronic phases in this material. In this study, we dis¬ 
covered that a moderate pressure (substantially below 
the critical pressure of the structural transition) not only 
closes the band gap, but also induces an electronic topo¬ 
logical transition (also known as Lifshitz transition |lClj ) 
that fundamentally modifies the material property of 
black phosphorus. Specifically, black phosphorus evolves 
into a semimetal above the Lifshitz transition. Conse¬ 
quently, we observed a colossal magnetoresistance up to 
approximately 80000% at a magnetic field of 9 T, which 
was possibly related to the nearly compensated electrons 
and holes in this semimetal state. Multiple electronic 


bands emerged from the Lifshitz transition, and produced 
quantum oscillations when the sample was subjected to 
external magnet fields. Analysis of the oscillation phases 
revealed that one of the emergent bands carried a non¬ 
trivial Berry’s phase, a prominent feature of a Dirac-like 
dispersion. 

We used single crystals of black phosphorus that were 
prepared using the high-pressure synthesis technique (de¬ 
tails provided in Supplemental Material [HI). The sam¬ 
ple was a p-type semiconductor at ambient pressure (Sup¬ 
plemental Material ll| ). The typical sizes of single crys¬ 
tal samples were 2.5x0.7x0.1 mm^, and the electrical 
contacts were arranged in a conventional Hall-bar config¬ 
uration. We applied hydrostatic pressure by using a self- 
clamped BeGu piston-cylinder cell that produces pres¬ 
sures up to 2.4 GPa. The magneto-transport measure¬ 
ments were performed in a Quantum Design PPMS-9. 

The application of hydrostatic pressure drastically 
changes the transport properties of black phosphorus, 
and induces a semiconductor-to-metal transition. Fig¬ 
ure 1 displays the temperature dependence of the lon¬ 
gitudinal resistivity p^x of a black phosphorus sample 
under varying pressures. The low-temperature resistive 
divergence was gradually suppressed by the applied pres¬ 
sure and the upturn feature completely vanished at P 
= 1.25 GPa. At P = 1.38 GPa, the high-temperature 
non-metallic behavior of pxx (T) was not detectable up 
to room temperature. Metallic resistivity in the entire 
temperature range from 2 to 300 K was observed at 
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higher pressures. Such behaviors can generally be cap¬ 
tured by the ratio between the resistivity at room temper¬ 
ature and that at the lowest temperature, also referred 
to as the residual resistivity ratio (RRR). As shown in 
the inset of Fig. 1, a prominent kink in the RRR indi¬ 
cates a semiconductor-to-metal transition occurring at a 
critical pressure of Pc — 1.2 GPa. Intriguingly, such a 
critical pressure coincides with a drastic enhancement of 
the sample magnetoresistance (MR) when the magnetic 
field is applied normal to the P atom plane (along the 
c-axis). The concurrent sharp enhancement of RRR and 
MR points to a profound change in the electronic struc¬ 
ture of black phosphorus at P = Pc- 



FIG. 1: (color online). Longitudinal resistivity p^x of black 
phosphorus versus T at various hydrostatic pressures. Num¬ 
bers beside the curves represent applied pressures in giga- 
pascals. Inset: The value of the residual resistivity ratio 
Pa;a:(300 YC)jpxx(^ K) and the transverse MR ([(pa:2;(9T)- 
pxx{H=Q)\/pxx{H=Q)) with if || c at T = 2 K as a function 
of pressure. 


Previous studies have confirmed that the orthorhom¬ 
bic (A17) structure of black phosphorus is preserved up 
to 4.5-5 GPa without any discontinuity in the struc¬ 
tural parameters [iH [^ . meaning that the transition 
we observed was an isostructural electronic transition. 
Moreover, the slope change of the pressure dependent 
resistivity at the critical pressure (Supplemental Mate¬ 
rial ll|) was consistent with an electronic topological 
transition The field dependence of the Hall re¬ 

sistivity provided strong evidence of co-existing electron 
and hole bands above Pc, as shown in Fig. 2(a), a hall¬ 
mark of multiband conduction, became apparent at P > 


1.20 GPa and extended to high field under higher pres¬ 
sure. Furthermore, the sign of p^y reversed from negative 
(electron-type) to positive (hole-type) at H ~ 0.11 T at 
1.20 GPa, and the sign reversal field shifted upward as 
the pressure increased (Fig. 2(a), inset). The Hall re¬ 
sistivity changes its sign in magnetic field when one of 
the two types of carriers (electrons and holes) has larger 
carrier population and the other shows higher mobility, 
as observed in graphite [i^ and NbSb 2 Therefore, 
at Pc ~ 1.20 GPa, black phosphorus is likely to enter a 
semimetal phase in which both electrons and holes are 
present at the Fermi level. Based on all of these findings, 
we identify the semiconductor-to-semimetal transition in 
black phosphorus as the topological Lifshitz transition 
of the band structure: The applied pressure pushed the 
conduction and valence bands toward each other, and 
eventually produced a semimetal state when the bands 
crossed the Fermi level successively. 

A notable consequence of this electronic transition is 
the observation of a colossal MR in the semimetallic 
phase of black phosphorus. When the pressure increased 
from Pc ~ 1.2 GPa up to approximately 2.0 GPa, the pos¬ 
itive MR was enhanced by a factor of ~ 40 and reached 
the exceedingly high value of ^80000% in a magnetic 
field of 9 T, with no indication of saturation [Fig. 2(b)]. 
A similar non-saturating large positive MR has been re¬ 
ported in a number of semimetals 28, and at¬ 


tributed to electron-hole compensation 
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Bl. The 

colossal MR was also observed with H parallel to the 
basal plane (Supplemental Material ll|) in spite of the 
expected anisotropy i El¬ 


la) 


(b) 




FIG. 2: (color online), (a) Hall resistivity pxy of a black 
phosphorus single crystal as a function of the magnetic field 
at various pressures at T = 2 K. Inset: Low-field sign reversal 
of pxy at P > 1.20 GPa. The Hall data were antisymmetrized 
with respect to positive and negative magnetic fields to re¬ 
move the pxx component, (b) Field dependence of transverse 
MR for P > 1.1 GPa with P || c, and current / applied in the 
basal plane, at T = 2 K. 


Quantum oscillations in magnetic field serve as a pow¬ 
erful tool for probing the Fermi surface (FS) topology 
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across the pressure-induced Lifshitz transition. Shub- 
nikov de Haas (SdH) oscillations emerged even before 
the low-temperature upturn of resistance was completely 
suppressed. Hints of SdH extrema on MR curves were al¬ 
ready discernible at P — 0.65 GPa (arrows in Supplemen¬ 
tal Material Fig. S2 llj) and the oscillation frequency 
could be resolved above 0.75 GPa. Oscillation patterns 
at various pressures are plotted in Fig. 3(a). Glear peaks 
manifested in fast Fourier transform (FFT) spectra [Fig. 
3(b)]. The oscillation frequencies are remarkably low, re¬ 
flecting the small sizes of Fermi pockets. At pressures 
below 1.25 GPa, a magnetic field of 9 T can force all the 
electrons into the lowest Landau level. According to the 
Lifshitz-Onsager quantization rule, the quantum oscilla¬ 
tion frequency (F) is proportional to the extremal cross- 
sectional area Sp of the FS [l^, F = {h/2TTe)SF . The 
lowest frequency measured with H || c was F = 4.0 T at P 
= 0.92 GPa, corresponding to a cross-sectional area Sp = 
0.038 nm“^, which was only about 0.014% of the area of 
Brillouin zone in (001) plane. More than one oscillation 
frequency could be resolved at F > 1.25 GPa (Fig. 3(b)), 
indicating the existence of multiple FS extrema. No¬ 
tably, the starting point of the multiple-frequency region 
was consistent with the onset pressure of the semimetal 
phase. The emergence of a new SdH branch strongly sug¬ 
gests that the band topology is modified at this electronic 
transition: Gombining this phenomenon with the Hall re¬ 
sult, there is a reasonable assumption that an electron FS 
appears at the critical pressure. 

In the FFT spectra of the semimetal phase, two sets 
of spectrally prominent features, a and /3, were identi¬ 
fied based on their pressure dependence above 1.38 GPa 
[Fig. 3(b) and (c)]. An Expansion of the FSs under 
pressure could be observed in both branches. At high 
pressures, branch a developed into the strongest feature, 
whereas branch /3 became comparably weak. We note 
that multiple peaks appear in the FFT spectra between 
1.25 GPa and 1.38 GPa, suggesting the existence of ad¬ 
ditional FS branches. The quantitative identification of 
such branches is however difficult in the data set currently 
available. 

Analyzing the temperature dependence of the oscilla¬ 
tion amplitude based on the Lifshitz-Kosevich (LK) for¬ 
mula (23| yielded small carrier cyclotron masses: the cy¬ 
clotron mass ratio m* = mjme (where m^ is the free 
electron mass) varied from 0.043 to 0.136 in the pressure 
range from 0.75 to 1.59 GPa and its pressure dependence 
mostly matched that of the corresponding oscillation fre¬ 
quencies [Fig. 3(d)] (for details see Supplemental Ma¬ 
terial [HI). The extracted carrier masses were consid¬ 
erably smaller than that reported for ambient pressure 
d, m, Q and this notable mass reduction reflected a 
drastically modified electronic structure under pressure. 
SdH oscillations were also clearly resolved when a mag¬ 
netic field was applied in the basal plane, and small Fermi 
pockets and small cyclotron masses were again observed 


(Supplemental Material Fig. S7 [ll|). The oscillations 
observed for F T c confirmed that at least one of the FSs 
was three dimensional (3D) in nature with close orbit 
perpendicular to the basal plane. 
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FIG. 3: (color online), (a) Negative second derivative of re¬ 
sistivity with respect to the magnetic field, -d^pxx/dlfi, at 
varying pressures. Magnetic field H was applied along the c 
direction. Data were collected for a black phosphorus single 
crystal labeled sample A. (b) The FFT spectra of -d?p^x/dlP 
of sample A. The labels a and /3 are assigned to the resolved 
FFT peaks, (c) The oscillation frequencies and correspond¬ 
ing cross-sectional areas of Fermi surfaces of sample A (solid 
symbols) and sample B (open symbols), plotted as a function 
of P. The uncertainty is determined by the half width at half¬ 
height (HWHH) of the FFT peaks. The raw data obtained 
for sample B is shown in Supplemental Material Fig. S6 [ll|. 
(d) Pressure dependence of the cyclotron mass ratio m* = 
m/me. m* was calculated based on the LK formula fitting of 
the T dependence of the normalized amplitude of the oscilla¬ 
tory component Ap^x/pxx (solid symbols) or the normalized 
FFT amplitude (open symbols). The average inverse fields in 
the FFT interval are noted beside data points. 


We now turn to the phases of quantum oscillations 
in the pressure-induced semimetal state. The SdH 
oscillation can be semi-classically described as follows 


F 

Aaxx oc cos[27r(— - 7 -b 5)], 


( 1 ) 


where a^x is the longitudinal conductivity, F the oscilla¬ 
tion frequency, 7 the Onsager phase factor and 6 an ad- 
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ditional phase shift that takes a value between - 1/8 and 
- 1 - 1/8 depending on the degree of three-dimensionality of 
the FS. The Onsager phase 7 relates to the Berry’s phase 
4>b through 7 = l/2-(/)B/27r [^, 38|. In an ordinary elec¬ 


tron system that exhibits parabolic energy dispersion, 
is zero and therefore 7 = 1/2. By contrast, in a mass¬ 
less Dirac electron system showing linear energy disper¬ 
sion, a non-trivial Berry’s phase = tt exists, yielding 
7 = 0 [ii, [a^. For the frequency component a which 
is spectrally dominate above 1.6 GPa [Fig. 4(a)], we 
could extract the value of 7 from the Landau index plots 
[ 3 ^ . The linear extrapolation of 1/versus the integer 
n plot is presented in Fig. 4(b), and the inset summa¬ 
rizes the - 7 -I -15 values for the a pocket, all of which are 
close to 0 with offsets lower than 0.07. Since black phos¬ 
phorus is a quasi-two-dimensional system, the absolute 
value of the phase shift 5 can be smaller than that in the 
3D limit (1/8) [i^. Thus the result of intercept anal¬ 
ysis indicated a phase factor 7 ~ 0. The zero Onsager 
phase corresponded to a non-trivial Berry’s phase 
TT, providing evidence of the linear band dispersion [ 2 ^ 
for pocket a (analysis of the intercept for other pock¬ 
ets is less conclusive owing to the relatively weak signal). 
Our results suggested that black phosphorus could be a 
new candidate for the elemental semimetals that possi¬ 
bly host exotic electronic states, besides graphite 4ll. l42| 
and bismuth [43j| . 


(a) 






FIG. 4: (color online). (a) Oscillation patterns of the 
magneto-conductivity axx{H) of pocket a, obtained by sub¬ 
tracting a smooth background from o-xx=Pxx/{Pxx+Pxy)j at 
various pressures with H \\ c. Integer indices are marked at 
the maxima of Aaxx- (b) Landau index n versus I/^qH for 
the SdH oscillation of pocket a. The linear extrapolations of 
the indices are indicated by solid lines. Inset: Intercepts of 
the linear extrapolations on n-axis, which equal to the total 
SdH phase -y-l-J. The correction terms <5 in the 3D limit, 5 = 
-1/8 (FS maxima) and 5 — -1-1/8 (FS minima), are denoted 
by dashed lines for the phase of Dirac electrons (7 = 0). Error 
bars represent the uncertainties of linear fitting. 


The Lifshitz transition occurring at Pc — 1.2 GPa 
stems from the lattice distortions of black phosphorus un¬ 
der pressure. Although the covalent bonds between phos¬ 
phorus atoms resist against compression, the change in 
the bond angle and the weak van der Waals interlayer in¬ 
teractions result in large compressibility in a and c direc¬ 
tion, respectively Q. Therefore, black phosphorus is ex¬ 
tremely soft and exhibits an extraordinarily pronounced 
pressure-induced band transformation for an inorganic 
semiconductor. Specifically, a recent first-principles cal¬ 
culation showed that the orbitals of adjacent atoms 
in distinct sub-lattices of the puckered lattice of black 
phosphorus gradually overlap under pressure (d^ j. This 
orbital overlap drives the band crossover near the Z point 
in Brillouin zone, develops multiple Dirac cones and pro¬ 
duces nearly compensated electron and hole pockets. In 
their band structure calculation there are three types of 
Fermi pockets in the semimetal phase and two of them 
are close to Dirac points [4^. This prediction is par¬ 
tially supported by our observations, although the third 
FS and a second Dirac-like band remain to be identified. 
Similar trend of band structure evolution in black phos- 
horus has been suggested in previous theoretical studies 


The qualitative consistence of calculation works and 
our experimental results strongly suggests a plausible 
scenario: black phosphorus evolves from a narrow-gap 
semiconductor to a 3D Dirac semimetal under external 
pressure, as a result of the formation of band crossover 
with linear dispersions. The Lifshitz transition we ob¬ 
served can be related to the theoretically predicted split¬ 
ting of an initial four-fold degenerate Dirac point at the 
Z point 4^. An alternative explanation is that the van 
Hove singularity of the upper band passes the Fermi level 
at the critical pressure [l9j |. resulting in a new electron 
pocket. We note that one FS emerges below Pc as an 
SdH frequency was already resolved at ^ 0.7 GPa. Since 
the black phosphorus samples are slightly p-doped and 
the Fermi level is closer to the top of valence band, it 
is possible that before the complete closure of the band 
gap, the Fermi level crosses the edge of valence band and 
a small hole pocket is formed. Therefore, when pressure 
is applied, black phosphorus may undergoes more than 
one electronic transition. Additional investigations are 
required for determining the exact pressure at which the 
band gap closes and for clarifying the details of band 
structure evolution under pressure. 

In summary, we observed a phase transition that in¬ 
troduces a typical semimetal phase in the direct gap 
semiconductor black phosphorus under a moderate hy¬ 
drostatic pressure of Pc — 1.2 GPa, which we identified 
as a topological Lifshitz transition of the electronic band 
structure. The semimetal phase was characterized by the 
colossal MR, the coexistence of electron and hole pock¬ 
ets, low carrier densities, and small cyclotron masses. 
A non-trivial Berry’s phase was detected for one of the 
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pockets, indicating a potentially peculiar band topology. 
Our results demonstrate that pressure, apart from layer 
number, serves as an effective parameter of tuning the 
electronic structure of black phosphorus. The dramatic 
response of black phosphorus to modest pressure reflects 
this materials capacity for sensitively detecting pressure 
or strain. Moreover, the evidence of linear dispersion in 
bulk black phosphorus under pressure points to a poten¬ 
tial two-dimensional Dirac fermion system in phospho- 
rene, a monolayer of black phosphorus. 
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